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Algorithms  for  solvents  and  spectral  factors  of 
matrix  polynomials 

LEANG  S.  SHIEHf,  YIH  T.  TSAYt  and  NORMAN  P.  COLEMAN  J 

A  generalized  Newton  method,  bated  on  the  contracted  gradient  of  a  matrix  poly¬ 
nomial,  is  derived  for  *olving  the  right  (left)  solvents  and  spectral  factors  of  matrix 
polynomials.  Two  methods  of  selecting  initial  estimates  for  rapid  convergence  of  the 
newly  developed  numerical  method  are  proposed.  Also,  new  algorithms  for  solving 
complete  sets  of  the  right  (left)  solvents  and  spectral  factors  without  directly  using  the 
eigenvalues  of  matrix  polynomials  are  derived.  The  proposed  computer-aided 
method  can  be  used  to  determine  the  spectral  factorization  of  a  matrix  polynomial  for 
optimal  control,  filtering  and  estimation  problems. 


1.  Introduction 

Linear  time-invariant  multi-input  multi-output  systems  are  often  described 
by  a  set  of  coupled  high-degree  multivariable  differential  equations.  The 
Laplace  transform  descriptions  of  such  systems  are  the  matrix  fraction 
descriptions  or  polynomial  matrix  descriptions  (Kailath  1980).  The  non¬ 
singular  denominator  matrix  polynomial  of  the  matrix  fraction  description, 
called  a  lambda  (A)  matrix  (Frazer  et  al.  1952,  Lancaster  1966),  characterizes  the 
properties  of  a  multivariable  (multi-inpuV-output)  system.  The  stability 
property  of  a  A-matrix  has  been  investigated  by  Lancaster  (1966),  Jury  (1974), 
Anderson  and  Bitmead  (1977),  Papaconstantinou  (1975),  Shieh  et  al.  (1978), 
Shieh  and  Sacheti  (1978),  and  Shieh  and  Tajvari  (1980).  Recently,  the  algebraic 
theory  of  A-matrices  has  been  investigated  by  Dennis  et  al.  (1976),  Gohberg  et  al. 
(1978),  Markus  and  Mereuca  (1978),  and  Tsay  et  al.  (1981).  Also,  numerical 
methods  for  so  lying  solvents  and  spectral  factors  of  a  matrix  polynomial  have 
been  proposed  by  Dennis  et  al.  (1978)  and  Shieh  and  Chahin  (1981). 

In  this  paper,  a  new  numerical  method,  based  on  the  contracted  gradient  of  a 
matrix  polynomial,  is  proposed  for  solving  a  complete  set  of  right  (left)  solvents 
and  spectral  factors  of  a  A-matrix.  The  proposed  method  can  be  applied  to 
determine  the  spectral  factorization  of  a  matrix  polynomial  for  optimal  control, 
filtering  and  estimation  problems. 

2.  Solvents  and  spectral  factorization  of  A-matrices 

An  nth-degree  mth -order  monic  A-matrix  is  defined  as 

A(A)=  £  At A»*<  (1) 

i-0 

where  Ait  t»l,  ...,n  are  mxm  complex  matrices  and  A0(,**lm)  is  an  m  x m 
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identity  matrix.  A  is  a  complex  variable.  Let  X  be  an  m  x  m  complex  matrix. 
Two  matrix  polynomials,  defined  by 

£  A{X-<  (2) 

t-0 

and 

Al(X)-  £  X~*At  (3) 

i-0 

are  referred  to  as  the  right  and  left  matrix  polynomials  associated  with  A-matrix 
.4(A),  respectively.  A  right  solvent  R  of  4(A)  is  defined  by 

4r(J?)«  £  AiR*-i~Om  (4) 

•  -0 

and  the  left  solvent  L  of  4(A)  is  defined  by 

4l(L)-  £  2>"<4(« 0.  (5) 

»-o 

where  0m  is  an  m  x  m  nail  matrix.  The  right  (left)  solvents  play  an  important 
role  in  the  analysis  of  the  A-matrix.  An  important  set  of  solvents  is  the  complete 
set  of  regular  solvents  (Gohberg  et  al.  1978,  Markus  and  Mereuca  1978).  The 
complete  regular  right  solvents  Rt,  i  —  1, ....  »  have  the  following  properties  : 

*Vi,  *,j«l, ....  n  (So) 

and 


u  p(J2«)-p(4( A))  (6  b) 

i-i 

where  p(Rj  is  the  spectrum  of  Rt  and  p(4(A))  is  the  spectrum  of  4(A). 

Similarly,  a  complete  set  of  regular  left  solvents  Lv  »«1, ...,  n  has  the 
following  properties : 

pfAJOp^)-^  *  ft  j,  ....  n  (7«) 

U  p(L<)*=p(4(A))  (7  6) 

i-l 


The  existence  and  uniqueness  of  the  complete  set  of  solvents  have  been 
investigated  by  Lancaster  (1986),  Dennis  et  al.  (1976),  Gohberg  et  al.  (1978) 
and  Markus  and  Mereuca  (1978). 

If  a  monic  A-matrix  can  be  decomposed  into  the  product  of  first-degree 
linear  A-matrices : 

A(X)^(XIn-Slt)(Xlm-Sn.t)...  (A 7.-5,)  (8o) 

then  (XJm-Si),  t»l, .... «,  are  referred  to  as  the  complete. set  of  linear  spectral 
factors.  The  mxm  complex  matrices  Sit  »  —  1, ...,  n,  are  the  spectral  factors 
of  the  A-matrix  4(A).  Notice  that  5,  is  a  right  solvent  of  4(A),  whereas 
is  the  left  solvent  of  4(A).  Other  spectral  factors  are,  in  general,  not  the  right 
or  left  solvents  of  4(A). 


or 
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Since  is  a  right  solvent  ( A  ifx)  of  /1(A),  a  deflated  A-matrix  can  be  obtained 
by  factoring  out  the  linear  spectral  factor  (XI  „  -  Rx)  from  A( A),  as  follows  : 

A<»(h)*A(X)(Xlm-B1)-*  =  (\lm-Sa) ...  (A Jm-St)  (8b) 

Determining  a  right  solvent  ( of  A  (1)(A)  yields  the  spectral  factor 
St  in  (8  b).  By  continuing  the  same  procedure,  we  can  determine  all  spectral 
factors  of  A(X)  in  (8  a).  Notice  that  Rta)  is  a  right  solvent  of  AU>(A)  but  not 
that  of  A(  A).  In  the  same  fashion,  the  spectral  factors  of  A(A)  can  be  computed 
from  the  successive  left  solvents  of  A(A)  and  the  deflated  .4(A). 

In  optimal  control  and  estimation  problems,  .4(A)  is  an  even  degree  real 
A-matrix  with  the  property 

A(A)«  AT(  —  A)  (9  a) 

where  T  designates  transpose.  The  spectral  factorization  (Youla  1961, 
Anderson  1967)  of  A  (A)  is  described  by 

AUHA^AM^-A)  (9  6) 

where  Ax(  A)  contains  all  latent  roots  of  .4(A)  in  the  LHP  of  the  complex  plane  and 
A,(  —  A)  in  the  RHP  of  the  complex  plane.  To  determine  'Aj(A)  the  procedure 
for  obtaining  the  spectral  factors  of  A(A)  in  (8a)  can  be  applied.  However,  by 
using  the  property  shown  in  (9  a),  the  computations  can  be  simplified.  In 
other  words,  if  R  is  a  right  solvent  of  A  (A)  in  (9)  with  all  eigenvalues  in  the 
RHP,  then 

A(A)  =  Aa>(A)(A/m  -  R)  (10a) 

where  A11 ’(A)  is  an  (»—  1)  degree  A-matrix.  From  the  property  of  (9  6)  we  have 

A( A)  »*  AT(  -  A) « ( -  XI m  -  RT)A  «>T(  -  A)  (10  6) 

From  (10  6)  we  observe  that  —R1  is  a  left  solvent  of  A  (A).  Thus,  the  factoriza¬ 
tion  of  A(A)  in  (9  6)  can  be  simplified  by  factoring  out  (A Im—R)  and  (XI  n  +  -RT) 
from  the  right  and  left  of  A(  A),  respectively,  if  a  right  solvent  R  is  found.  The 
deflated  A(A)  becomes 

A»\X)m(XIm+Rr)-iA(X)(XIm-R)-1  (10  c) 

where  A0)(A)  is  an  (n  -  2)  degree  A-matrix. 

Repeating  the  same  procedure  for  (10  c)  gives  the  desired  spectral  factors. 
Thus,  the  spectra]  factorization  of  A(A)  can  be  determined  from  the  obtained 
spectral  factors. 

In  this  paper,  a  new  numerical  method  for  solving  the  right  (left)  solvents 
Ri(Lt)  and  spectra]  factors  (S{)  of  A  ( A)  is  derived  by  using  the  contracted  gradient 
of  a  matrix  polynomial. 

3.  Gradient  of  a  matrix  polynomial 

The  gradient  of  an  nth-degree  mth-order  matrix  polynomial  is  defined  as  a 
fourth-order  tensor  VA(X)  with  components  (Vetter  1973,  Fong  1971) 

{VA(J)}<./.*.f-~-{A(Z)}u;  i,j,k,l-l,...,m  (11) 

3c  2 
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A 


where  Xkt ,  denotes  the  (i,  I)  element  of  X ,  {A  (X)}u  designates  the  (i,  j)  element 
of  and  {VJ4(Z)}<iy  iil  denotes  the  (t,  j,  k,  l)  element  of  VA(.X). 

Direct  use  of  the  gradient  for  the  purposes  of  this  paper  involves  the  inversion 
of  a  fourth-order  tensor,  thus  causing  computational  difficulty.  A  contraction 
operation  (Fong  1971)  on  VA(.X)  with  respect  to  an  arbitrary  mxm  square 
matrix  Y  is  defined  as 

I  I  *’.i- i . *  d*«) 

i-J  1-1 

from  which  it  follows  that 

(12fc) 

where  ij  is  a  scalar  variable.  Having  the  contraction  operation  defined  in  (12), 
the  contracted  gradient  V  A(X)Y  becomes  an  mxm  matrix.  Thus,  compu¬ 
tational  algorithms  can  be  easily  derived  and  the  computational  difficulty  can  be 
greatly  reduced.  The  contraction  operation  in  (12)  can  also  be  extended  to  the 
right  (left)  matrix  polynomial  in  (2)  and  (3),  as  follows. 

Given  a  A-matrix  in  ( 1 )  and  the  right  matrix  polynomial  in  (2),  the  contraction 
operation  on  VAR(Z)  with  respect  to  a  matrix  Y  is 

vAB(d)y-|.(AR(i+i,r))U 

-  f  1  (i.d+^nu  (i3a> 

i-Q 

Each  term  inside  the  summation  of  (13  a)  can  be  computed  as 

1  (At(X +,DT)U-^«  l‘  X'YX'-*-'  (13  6) 

d7f  f-0 

Substituting  (13  b)  into  (13  a)  and  rearranging  the  indexes  gives 

var(Z)f-  "£  (i*) 

imO  f — 0  * 

Performing  index  transformations  or  letting  k-i+q+ 1  and  j  * i,  ( 1 4)  becomes 
VAR(J)y»  f  BkK&)YX'~*  (16  a) 

t-i 

where  BkR(X)  is  the  right  matrix  polynomial  of  the  following  A-matrix  Bk(  A) : 

S*(A)«  *1  A,**-’-1  (16  6) 

i-o 

Similarly,  for  the  left  matrix  polynomial  AL(X),  we  have 

VAL(Z)y.  i  X"*YBkh(X)  (16) 

t-i 

where  Bkh(X)  is  the  left  matrix  polynomial  of  the  A-matrix  Bk(  A)  defined 
in  (18  b). 


1 
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4.  A  generalized  Newtoa  method  for  solving  solvents 

The  right  matrix  polynomial  4H(X)  can  be  expanded  around  an  m  x  m 
matrix  Xt  as 

AR(X)-AR(X0)  +  VAR(XO)(X  -X,)  +  0((J i-X0))  (17) 

where  O((X~X0))  is  a  matrix  polynomial  with  high-degree  terms  of  (X  -  X0), 
and  V.4R(X0){X  -  X0)  is  a  contracted  gradient  of  dimension  »t  x  m. 

The  first-degree  approximation  of  (17)  with  |  AXjj  <  1  becomes 

AR(X0)  +  VAR(X0)AX  (18) 

where  AX  AX-  X9. 

Define  a  recursive  formula  * 

•^t+i  “  X<  +  (19a) 

Thus,  (18)  becomes 

(19  6) 

If  XM  is  the  right  solvent  of  .4(A),  or  -4R(Xi+l)«0m,  then 

^(XJ  +  V^XJAX^.O.,  (20) 

Solving  AX1+1  from  (20)  and  substituting  it  into  (19  a)  gives  the  recursive  formula 
for  solving  the  right  solvent  of  4(A).  To  solve  the  AX<+1  we  use  the  contracted 
gradient  derived  in  (18),  as  follows. 

From  (15)  we  have 

VAR(XJAXi+1  -«  ^kRiX^AXi+iXi"-*  (21) 

i>i 

Substituting  (21)  into  (20)  yields 

Z  ^kR(-^i)AXi+lXin~lc=  —  AR(Xt)  (22) 

*-i 

AXi+1  can  be  solved  by  using  Lancaster’s  Kronecker-product  expansion  method 
(Tsay  et  al.  1981,  Lancaster  1970)  as  follows  : 

Vec  {AX<+1}  -  -  ^(X,)-1  Vec  MR(Xf)}  (23  a)  ' 

where 

<W«  f  (^--*)t®^r(^)  (23  ft) 

®  designates  the  Kronecker  product  (Barnett  1971)  and  Vec  {Jl/}  is  a  column 
vector  consisting  of  the  column  vectors  M}  of  an  m  x  m  matrix  M  or 

Vec{Af}AlAf1T,JItT,...,MmT]T  (23  c) 

Note  that  the  dimension  of  Q(X{)  is  m*  x  m*.  BkR(X{)  is  a  constant  matrix  of 
the  right  matrix  polynomial  of  the  A-matrix  2?*(A)  shown  in  (15  ft)  evaluated 
at 

Similarly,  the  recursive  formula  for  solving  the  left  solvents  of  A( A)  is 


(24) 
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where  AXi+1  is  the  solution  of  the  following  linear  matrix  equation  : 

V4l(Z,)AX,+1-  ^  BmAXJs*  —  A-^Xi)  (25) 

*-i 

or 

Vec  {AJ<+1}  -  -  Vec  {4L(Z{)}  (26  a) 

where 

*(*,)«  2  (26  5) 

The  convergent  criterion  in  both  (19)  and  (24)  is 

\AXm\«  (27) 

where  « is  an  assigned  small  value. 

Note  that  when  m  —  1,  .4(A)  in  (1)  becomes  a  scalar  A- polynomial  The 
recursive  algorithm  in  (19)  and  (24)  is  the  well-known  Newton  method.  Thus, 
the  proposed  algorithm  in  this  section  is  a  generalized  Newton  method. 

5.  Algorithms  for  solving  complete  sets  of  sotants 

To  construct  a  similarity  block  transformation  (Dennis  et  al.  1976,  Tsay  et  al. 
1981)  for  the  analysis  of  a  large-scale  multivariable  system,  it  is  necessary  to 
determine  a  complete  set  of  solvents. 

Given  a  A-matrix  4(A)  in  (1),  if  aright  solvent  if  is  obtained,  then  4(A)  can 
be  described  as 


4(A)  — 2J(A)(Af„— ilj) 

(28  a) 

where 

B( A)-  "£  -BiA-1-4 

(28  5) 

and  the  coefficients  B{  can  be  determined  by 

B% m 4{ + i  —  1,  • a— -1 

(28  c) 

The  right  solvent  of  B( A)  is  a  spectral  factor  of  4(A)  but  in  general  not  the 
right  solvent  of  4(A).  To  resolve  this  problem,  we  have  to  determine  the  left 
solvent  of  4(A)  associated  with  J?x  by  using  the  following  relationship  (Tsay 
H  at.  1981) : 

Lx  —  Q{~1  RiQx  if  rank  Qx 

(29) 

where  Q,  is  the  solution  of  the  following  linear  matrix  equation  (Tsay  et  al.  1981, 
Lancaster  1970) : 

2 
im  0 

(30) 

or 

Vec  {$*}-[  "f  BTgKV^^Vec  {/„} 

(31) 

Algorithms  for  matrix  polynomials 


1309 


Equations  (29)  and  (30)  can  be  easily  verified  by  substituting  £,  at  (29)  into  the 
left  matrix  polynomial  of  A( A)  in  (28),  or  Al(£i),  and  using  the  result  of  (30). 
Also,  following  the  similar  formulation  in  (25)  and  solution  in  (20),  we  can  solve 
(30)  for  Qv  which  is  shown  in  (31 ).  Once  is  found,  A{ A)  can  be  expressed  by 


A(  A)  -  (A/.-LJA^A)  -  AA  a)(A)  -i^A  *>(A) 

(32  a) 

A<»>(A)-  J  At«>A— ^ 

(32  b) 

<■0 

A^-Aj+^A*.^,  t-1,2, 

(32  c) 

Ammming  that  jR*  is  a  right  solvent  of  A0)(A),  or  A  a,(ff1)»0ml  then  from  (32  a) 
we  have  the  right  matrix  polynomial 

i8(121)-4H«i(il)Jf, -InAj.0  >(*,)  (33) 

Therefore,  a  right  solvent  of  A0)(A)  is  a  right  solvent  of  A  (A).  Repeating  the 
same  procedures,  we  can  determine  the  complete  set  of  right  solvents.  The 
algorithm  can  be  summarized  as  follows. 

Algorithm  1 — Complete  set  of  right  solvents 

(1)  Set  an  index  i  —  1. 

(2)  Find  a  right  solvent  Ht  of  A(X)  by  using  the  recursive  algorithm  in  (19  o) 
and  (23). 

(3)  Find 

5(A)- 

where 

Bt  -  A,  +  B}^iRu  j  —  1, 2, ....  *— *,  and  5#— A, 

(4)  Find  L(  -  G<-1  •*<{?<  ^7  u«ng  (31) 

Vec  {<?«}- j  Y  5/®(5<-<-o}“1  Vec  {/„} 

(5)  Construct 

Aa,(A)-  £  A,  A—*-* 

where 

j  —  1,  2,  ....  «  — i,  and  A0G>-A, 

(6)  Replace  A(A)  by  Aa,(A)  and  raise  the  index  t  by  1. 

(7)  If  i< n,  then  go  to  (2) ;  otherwise  stop. 

Similarly,  the  algorithm  for  determining  the  complete  set  of  left  solvents  is 
described  in  Algorithm  2  as  follows. 


Algorithm  2— Complete  set  of  left  solvents 

(1)  Set  an  index  i- 1. 

(2)  Find  a  left  solvent  Li  of  .4(A)  by  the  recursive  algorithm  in  (24)  and  (20). 

(3)  Find 

5(A)-  Y 

J- o 

where 

Sj^Aj- )  *1| ...» s  —  and  £|s4| 

(4)  Find 

Bt-PtLiPc* 

where  the  mxm  matrix  P4  can  be  constructed  by 


Vec 


{5<}«{  *f  (i^)T®^}_1Vec  {/„} 


(5)  Construct  -da,(A)  of  degree  »-*  by 


»-« 


A«>( A)-  £  Aia) 

i~0 


where 


+  j-1,  and  A,«>-A0 

(6)  Replace  .4(A)  by  A(1)(A)  and  raise  the  index  i  by  1. 

(7)  If  t  ^  n,  then  go  to  (2) ;  otherwise  stop. 

If  a  set  of  linear  spectral  factors  is  desired,  then  the  following  algorithm  which 
is  derived  from  (8  a)  and  Algorithm  1  can  be  applied. 


Algorithm  3 — Linear  spectral  factors 

(1)  Set  an  index  i«  1. 

(2)  Find  a  right  solvent  5<(  *  S{)  of  A  ( A)  by  the  recursive  algorithm  in  (19  a) 
and  (23). 

(3)  Find 

5(A)-  Y  5^— *-* 

J-o 

where 

J5y- Aj  +  Bj.iSi,  j-1, ....  n  —  i,  and  2?S-At 

(4)  Replace  A  (A)  by  R(A)  and  raise  the  index  f  by  1. 

(5)  If  i  <  «,  then  go  to  (2) ;  otherwise  stop. 

When  the  A-matrix  of  interest  has  the  property  that  A(A)-AT(  — A)  of 
2n  degree  and  the  spectra]  factorization  is  desired,  we  have  the  following 
algorithm. 
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Algorithm  4 — Spectral  factorisation 

(1)  Set  an  index  i  >■  l. 

(2)  Find  a  right  solvent  R{{  —  S{)  of  A  ( A)  by  the  recursive  algorithm  in  ( 1 9  o) 
and  (23). 

(3)  Find 

n-tf  +  l 

B(  A)«  £  »  a— M+w 

J-O 

where 


Bj ■^1^4'  j *  1| ....  n  — 2»+ 1,  and  ^ 

(4)  Find 

where 

Af'-Bj-SJ  Bs_lt  j«l, .... »  — 2»,  and  4,a>-210 

(5)  Replace  4(A)  by  4a>(A)  and  raise  the  index  t  by  1. 

(6)  If  i  4  n/2,  then  go  to  (2) ;  otherwise  stop. 

The  resulting  spectral  factorization  is 

-4(A)  -  Lft(  -  A)-D(A) 

where 

D(X).(Xlm-SK/t) ...  (XIm-8t){\Im—S1) 


6.  Normalization  and  initial  guesses 

The  rate  of  convergence  of  the  proposed  numerical  method  heavily  depends 
upon  the  initial  guesses.  In  order  to  derive  a  systematic  method  for  deter¬ 
mining  the  initial  guesses,  the  A-matrix  is  normalized  such  that  the  latent  roots 
of  the  A-matrix  are  uniformly  distributed  around  a  unit  circle  with  centre  at 
A-0. 

The  arithmetic  mean  of  the  latent  roots  of  a  given  A-matrix  in  (1)  is 


£  A<  —  trAl 

i-i  _ 

nm  n 


(34) 


where  A*  are  the  latent  roots  of  4(A)  -  0  and  trA1  is  the  trace  of  a  matrix  Av 
Substituting 

A-y  +  i,  (36) 

into  (l)-results  in 

^(1>(y)-^’*+^ifl,y"-,+  ...  +  4.<*>  (36) 

where 
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The  geometric  mean  of  the  latent  roots  of  A  n)(y)  is 

K -  (  jjj  -  (det  i.®')1'"  (37) 

where  at  are  the  latent  roots  of  An)(y)  —  0. 

Assuming  that  £,5*0  and  substituting 

(38) 

into  (36),  we  have 

A  a\z)  =  7mz*  +  A,®*  z*-1  +  ...  +A^.1«>z  +  A,«>  (39) 

where  |det  A,'*]  *•  l. 

A  «>(z)  is  referred  to  as  the  normalized  A- matrix  of  A  ( A).  The  initial  guesses 
for  Algorithms  1  to  4  can  be  set  in  the  following  two  ways. 

(1)  Lin's  method  (Lin  1943) 

The  scalar  version  of  Lin’s  method  is  extended  to  the  matrix  case  for  deter¬ 
mining  the  initial  guess.  The  initial  guess  for  the  recursive  algorithm  of  the  right 
solvent  of  A(t)(z)  in  (39)  is 

-Ro  =  {^— if  rank  (40) 

The  initial  guess  for  the  recursive  algorithm  of  the  left  solvent  of  A(,)(z) 
becomes 

<*>}-»  ifrankA^^-m  (41) 

(2)  Unitary  matrix  method 

Since  the  latent  roots  of  the  normalized  A-matrix  has  its  arithmetic  mean  at 
zero  and  the  geometric  mean  at  unity,  the  initial  guess  can  be  Bet  to  a  matrix 
with  the  magnitudes  of  all  eigenvalues  equal  to  unity.  The  matrix  is  a  unitary 
matrix. 

The  desired  solvents  or  spectral  factors  of  the  A-matrix  can  be  determined 
by  denormalizing  the  solvents  or  spectral  factors  obtained  from  A  B)(z)  via  the 
transformations  in  (38)  and  (36). 


7.  Examples 

Two  examples  are  illustrated  for  finding  the  solvents,  spectral  factors  and 
spectral  factorization  of  A-matrices.  The  computations  were  carried  out  on  a 
PDP-11/70  using  single  precision  floating  point  only. 

Example  1 

Given  the  A-matrix 

[4  2]  [  12  II]  f!9  14] 

A*+  I  A+ 1  (42) 

-2  7j  [-2  28 J  L16  36J 

find  the  complete  set  of  regular  right  solvents  and  its  linear  spectral  factors. 


Ift 
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Using  the  normalization  scheme  in  §  6,  we  have 

=  —1*83333 
4,-1-071383 

and  the  normalized  matrix  polynomial  Aa\z)  in  (39)  becomes 


AW(z)-Itzi  + 


1-40009  1-8666781  [6- 

z*  + 

-1-86678  1-40009  J  (_4-< 


46154  3-19447" 


64650  10-8176  J 
‘  3-48240  0-4517701 


•5263  1-65273 


The  initial  guess  by  using  Lin’s  method  is 


-0-0734806  0-00712910" 


•941503  -  0-155843 


Applying  Algorithm  1  to  (44)  and  using  the  initial  guess  J?0  gives  the  right 
solvent  of  Aa)(z)  at  the  fourth  iteration  with  an  error  tolerance  of  *- 10~*  in 
(27).  Denormalizing  the  obtained  solvent  by  using  (35)  and  (38)  yields  the 
desired  right  solvent 

[- 2-00000  -  0-37887  x  10~H 

(45  6) 

-1-00000  -  2-00000  J 

The  corresponding  left  solvent,  which  can  be  determined  from  (29)  and  (31),  is 

f  -2-46154  0-1153851 


1-84615  —  0-153846 J 

Now  factoring  out  (A It-Lx)  from  the  left  of  A  (A),  or  using  (32),  we  have 


•53846  2-11538 


AO>(A)-/tA*+ 


•84615  5-46154 


f|  f  7- 
A+ 

J  L10, 


76923  6-42308] 


(46  a) 


10-7692  15-6923 


The  normalized  A-raatrix  with  —  — 1-750  and  4,—  2-70594  is 


-0-724901  0-781756 


A<»(2)-If*  + 


-1-42138  0-72490 


1  ri-11162  0-371635" 

z+  (46  6) 

1*06632  1-25607 


The  initial  guess  by  Lin’s  method  is 


0-0474374 


1-21660  1 
0-652733J 
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The  right  solvent  of  Aa  >(2)  can  be  obtained  at  the  thirteenth  iteration  with  the 
same  error  tolerance  of  t  =*  10~*.  After  denormalization  we  have  the  desired 
right  solvent  of  AU,(A)  or  another  right  solvent  of  A  (A),  as  follows  : 

[-0-999976  1-49996' 

(46  d) 

-2-00004  —  2-00001  _ 

The  corresponding  left  solvent  is 

[-6-94525  3-35920* 

(46  e) 

-  9-64537  3-94526_ 

After  factoring  out  (\It  —  Lj  from  the  left  of  A  (1)(  A)  we  have 

[-  5-40679  5-47459' 

(47  a) 

- 13-4915  9-40680 

From  (47  a)  we  have  the  desired  right  solvent  of  A  '*>(  A)  or  the  last  right  solvent 
of  A(A)  aa 

T  5-40679  -  5-47459' 

(47  6) 

(_  13-491 50  -  9-40680 

The  complete  set  of  linear  spectral  factors  of  A  (A)  in  (42)  can  be  determined  by 
using  Algorithm  3.  The  result  is 

A(X)^(XIt-S9)(Xlt-S9n\It-S1)  (48) 

where 


Si 

S, 


I 

! 

[ 


-2-00000 

-  0-37887  x  lO-H 

-1-00000 

—  2-00000  J/ler*Uoo  number -4 

-0-99996 

2-9999n 

-1-00002 

—  2-00004 J  ll<nlU<m  number- 13 

-1-00004 

-  4-99991 1 

4-00002 

—  2-99986 J  luntjog  number  -  0 

It  might  be  interesting  to  investigate  the  relationship  between  the  eigenvalues 
of  A  ( A)  and  the  corresponding  solvents  and  spectral  factors.  The  eigenvalues  of 
A(A)  or  solvents  and  spectral  factors  are 


<*Sl)-v(Rl)--2,  -2 

«7(S,)-<r(Jft)--l-5±jv/ll  («) 

<t(5,)  =  <r(J?,)-  -2±jV7 
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From  (49)  we  observe  that,  although  we  have  repeated  latent  roots  with  a 
Jordan  chain  of  length  2  and  two  pairs  of  complex  roots,  the  corresponding 
solvents  and  spectral  factors  are  all  real.  This  peculiarity  has  important 
applications  in  the  analysis  and  synthesis  of  large-scale  multivariable  control 
systems. 


Example  2 

Given  a  A-matrix 

[0  —  6")  f-37  31 

As  +  A* 

6  Oj  [  3  -  i8j 

[  0  49]  I” 68  2" 

A+  (60  a) 

-49  OJ  [2  26 

and 

A(A)  =  AT(-A)  (50  6) 

it  is  desirable  to  find  the  spectral  factorization  of  A(X)  such  that 

4(A)-4,(AMit(-A)  (51) 

where  j4x(A)  contains  all  latent  roots  of  A  ( A)  in  the  LHP  of  the  complex  plane  and 
-  A)  in  the  RHP  of  the  complex  plane. 

Using  Algorithm  4  and  the  initial  guess  lt  (a  unitary  matrix),  we  have 

^(AMA/.-^K  A/,-5,)  (52) 

where 


[-0-99999  2-00001 

-1-00000  -  0-99999 

-3-99998  1-99996 

-0-99999  -2-99999 

Thus,  At(A)  becomes 

[4-99997  -3-99997]  T 1-99996  -7-99996' 

I  A+  I  (53) 

1-99999  3-99998J  [4-99996  0-99990 

Observe  that,  although  the  latent  roots  of  A  (A)  are  all  complex,  i.e. 

A--l±jV2;  $(-7  ±jV7);  i±iV2;  1(7±JV7) 

Sx  and  St  are  real  matrices. 


(54) 
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A  generalized  Newton  method,  which  is  based  on  the  contracted  gradient 
of  a  matrix  polynomial,  has  been  derived  for  solving  the  complete  sets  of  solvents 
and  spectral  factors  of  a  matrix  polynomial  without  directly  using  the  eigenvalues 
of  the  matrix  polynomial.  Lin’s  method  and  the  unitary  matrix  method  have 
been  proposed  for  determining  the  initial  estimates  of  the  solvents  and  spectral 
factors.  Thus,  the  convergence  speed  of  the  new  numerical  method  is  improved. 
The  proposed  method  can  be  used  to  determine  the  spectral  factorization  of  a 
matrix  polynomial  for  optimal  control,  filtering  and  estimation  problems.  The 
obtained  solvents  can  be  used  to  construct  a  block  similarity  transformation 
xngtnx  for  the  decoupling  of  a  large-scale  multivariable  control  system. 
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